arXiv: 1503.08162v3 [hep-ph] 6 Dec 2016 


The B — L Supersymmetric Standard Model 
with Inverse Seesaw at the Large Hadron Collider 

S. Khalil 1 and S. Moretti 2 

1 Center for Fundamental Physics, Zewail City for Science and Technology, 6 October City, 

Cairo, Egypt. 

3 School of Physics & Astronomy, University of Southampton, Highfield, Southampton, UK. 


Abstract 

We review the TeV scale B — L extension of the Minimal Supersymmetric Standard 
Model (BLSSM) where an Inverse Seesaw mechanism of light neutrino mass generation 
is naturally implemented and concentrate on its hallmark manifestations at the Large 
Hadron Collider (LHC). 


1 Introduction 

The solid experimental evidence for neutrino oscillations, pointing towards non-vanishing 
neutrino masses, is one of the few firm hints for physics Beyond the Standard Model 
(BSM). Neutrinos are strictly massless in the SM due to two main reasons: (*) the absence 
of right-handed neutrinos; (■ ii ) the SM has an exact global Baryon minus Lepton (B — L ) 
number conservation. However, the minimal extension of the SM, based on the gauge 
group SU (3)c x SU(2)lxU(1)y x U (1) b-l, can account for light neutrino masses through 
either a Type-I Seesaw or an Inverse Seesaw (IS) mechanism DG0 

In the type-I Seesaw mechanism right-handed neutrinos acquire Majorana masses at 
the B — L symmetry breaking scale, which can be related to the Supersymmetry (SUSY) 
breaking scale, i.e., 0(1) TeV, therefore the Yukawa neutrino coupling must be < 0( 10 -6 ) 
[3]. In contrast, in the IS case, these Majorana masses are not allowed by the B — L gauge 
symmetry and another pair of SM gauge singlet fermions with tiny masses (0(1) keV) 
must be introduced. In this case, there is no severe constraint imposed on the neutrino 
Yukawa couplings and the possibility of testing this type of model at the LHC is quite 
feasible. Moreover, in the limit of the above mentioned tiny mass, we have essentially 
massless light neutrinos. Therefore, such a small scale can be considered as a slight 
breaking of a global symmetry, hence, according to’t Hooft criteria, the smallness of this 

1 For old attempts at analysing a high scale B — L extension of the SM, see Ref. [ 31 . 
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scale becomes natural. One of these two singlet fermions couple to right-handed neutrinos 
and is involved in generating the light neutrino masses. The other singlet (which is usually 
called inert or sterile neutrino) is completely decoupled and interacts only through the 
B — L gauge boson, therefore it may account for warm Dark Matter (DM) [5j, see also 
Refs. [6l E]. 

In both scenarios, this B — L model induces several testable signals at the LHC involv¬ 
ing the new predicted particles: a Z' (neutral gauge boson associated with the U(1)b-l 
group), an extra Higgs state (an additional singlet state is introduced to break the gauge 
group U{1)b-l spontaneously) and three (Type-I) or six (IS) heavy neutrinos, zy,. (that 
are required to cancel the associated anomaly and are necessary for the consistency of the 
model). This is the setup for the non-SUSY sector of the B — L scenario, which is well 
established in the literature (see Refs. for a review of its main phenomenological 

manifestations). 

It is the purpose of this paper to review its Supersymmetric version, the BLSSM, 
particularly in the IS framework. The reason to concentrate on the IS version of the 
BLSSM is twofold. On the one hand, its phenomenological exploitation is more recent 
and less evolved with respect to the type-I Seesaw mechanism case, hence the need of 
reviewing its status in view of the ongoing Run 2 stage of the LHC. Secondly, the former 
appears to have some intriguing and testable features in the Higgs sector that the latter 
has not (as we shall try to emphasise), again, calling for a timely assessment of the LHC 
sensitivity to it, given the high priority that Higgs physics has for Run 2. 

The plan of the paper is as follows. In the next section we proceed to the construction of 
the BLSSM Lagrangian. In Sect. [3] we describe how dynamical Electro-Weak Symmetry 
Breaking (EWSB) occurs in the BLSSM whereas in Sect. [4] we introduce its particle 
spectrum. In Sect. [5] we study in particular the Higgs masses. Then, in Sect. [6j we 
describe the main manifestations of the BLSSM at the LHC. We conclude in Sect. 0 

2 Constructing the BLSSM 

The particle content of this model includes the following superfields in addition to those 
of the Minimal Supersymmetric Standard Model (MSSM): ( i ) two SM singlet chiral Higgs 
superfields Xi , 2 , whose Vacuum Expectation Values (VEVs) of their scalar components 
spontaneously break the U{1)b~l and is required to cancel the U(1)b~l anomaly; 
(ii) three sets of SM singlet chiral superfields, 1 — 1,2,3), to implement the IS 
mechanism (also in order to cancel the B — L anomaly). Tab. [I] provides the particle 
content of the SUSY version of the B — L model with IS (henceforth BLSSM-IS for short) 
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Superfield 

Spin 0 

Spin 4 

Generations 

SU(3) C x SU(2) l x U (l)y x U(1) B . L 

Q 

Q 

Q 

3 

(3.2, 1,1) 

d c 

d c 

d c 

3 

(3.1. 5 .- 5 ) 

u c 

u c 

u c 

3 

(1.3-1 ,-|) 

i 

l 

l 

3 

to 

1 

tO|h^ 

1 

I— 1 

e c 

e c 

e c 

3 

(l,l, l,l) 

l) C 

v c 

v c 

3 

(1,1, 0,1) 

si 

Si 

Si 

3 

(1,1, 0,2) 

S2 

s 2 

S 2 

3 

(1,1, 0,-2) 

H d 

H d 

H d 

1 

(1.2. -1,0) 

H u 

H u 

H u 

1 

(1,2.1,0) 

Xi 

Xi 

Xi 

1 

(1,1,0,1) 

X2 

X2 

X2 

1 

(1,1, 0,-1) 


Tabic 1: Chiral superfields of the BLSSM-IS and their quantum numbers under SU(3)c x 
SU(2) l x U(1) Y x U(1) B -l ■ 


as well as the different charge assignments of each Superfield. The superpotential of the 
leptonic sector in this model is given by m 


W =Y u tiqH u -Y d dqH d - Y e elH d + Y„ u 1 //„ + Y,vx 1 h + l*H u H d - ft'X .2 (1) 


In order to prevent a possible large mass term Msis 2 , we assume that the superfields 
u, xi ,2 and s 2 are even under matter parity, while .Si is an odd particle. By assuming 
a minimal Supergravity (mSugra) inspired universality of parameters at the scale of a 
Grand Unification Theory (GUT), we obtain that the SUSY soft breaking Lagrangian is 
given by 


£ soft = ml 

+ |X2 


+ 


\q\" + |«| + \d\" + \l\ + |e c | + \v c \~ + |«5i| + | £ 2 1 + \Hd\ + \H U \ + |xi| 


Y u A qH u u c + Y d A qH d d c + Y A lH d e c + Y A IH U D C + Y s a v c X yS 2 

1 


+ [BiiiHiHi + n' x 1 X 2 ) + h.c.) + -M 1/2 [g a g a + W a W a + BB + B'B' + h.c. 

where the trilinear terms are defined as (Y^ 4 )^ = ( YfA)ij. 

The B — L symmetry is radiatively broken by the non-vanishing VEVs (Reyf) = 4= 
(i = 1,2) while the EW one by the non-zero VEVs (Rei/.J) d ) = v Utd /\/2, with v = 
\J v u + v 2 — 246 GeV, v' = \Jvf + v 2 — 0(1) TeV and the ratio of these VEVs are 
defined as tan/3 = v u /v d and tan/3' = v' x jv' 2 [!). After B — L and EW symmetry breaking, 
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the neutrino Yukawa interaction terms lead to the following expression: 

C m = m D V L v c + M r v c S 2 + h.c., (2) 

where mo = ~^Y u v u and Mr = ^YgV^. In this framework, the light neutrino masses are 
related to a small mass term n s S 2 in the Lagrangian, with /j, s ~ 0( 1) KeV, which can be 
generated at the B — L scale through a non-renormalisable higher order term where 
M is the mass of a heavy state whose loop(s) or tree-level tadpole diagrams generate 
the corresponding higher order term (M ~ 10 6 in our case).Therefore, one finds that the 
neutrinos mix with the fermionic singlet fields to build up the following 9x9 mass matrix, 
in the basis ( v L , u c , S 2 ): 


( 

0 

rn D 

0 

\ 


T 

m D 

0 

Mr 


V 

0 

M T R 

l^s 

) 


(3) 


The diagonalisation of the mass matrix, Eq. ([3]), leads to the following light and heavy 
neutrino masses, respectively: 


m Vl = itidMr 1 ix s (Mr) 1 m/,, (4) 

m Vh = m UH , = sjM 2 R + m 2 D . (5) 

Thus, one finds that the light neutrino masses can be of order eV, with a TeV scale 
Mr , if fi s <C Mr, and a order one Yukawa coupling Y v . Such a large coupling is crucial 
for testing the BLSSM-IS and probing the heavy neutrinos at the LHC. As shown in 
[12] , the mixings between light and heavy neutrinos are of order (9(0.01). Therefore, the 
decay widths of these heavy neutrinos into SM fermions are sufficiently large. It is worth 
mentioning that the second SM singlet fermion, Si, remains light with mass given by 


m Sl =Hs- Oil) keV, (6) 

where Si is a sort of inert/sterile neutrino that has no mixing with the active neutrinos. 
It can therefore be a good candidate for warm DM as emphasised in Ref. [5]. 


3 Radiative B — L symmetry breaking in the BLSSM-IS 

The breaking of B — L can spontaneously occur through the VEV of the scalar held 
Xi ,2 or the right-handed sneutrino [1], depending on the initial values of the Yukawa 
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couplings and soft terms involved in the Renormalisation Group Equations (RGEs) of the 
parameters in the scalar potential V(xi,X 2 ), where 

V(XhX 2 ) = lA\Xi\ 2 + ^l\X 2 \ 2 ~ fJ>l(xiX 2 + h.c.) + -g BL (\xz\ 2 ~ \xi\ 2 ) j (7) 

where g\ 2 — m \ 12 + |A* 7 ] 2 , (J% = —B'g! and g B L is the gauge coupling of U(1) B -l- The 
stablitity condition of V{x\, X 2 ) is given by 


q 2 / 2 1 2 

2 /i 3 < /i] + /i 2 - 


( 8 ) 


The minimisation of this potential, = 0, 2 = 1,2, implies that 


M 2 


g\ = gl cot j3' H-cos 2/3', 


4 

M 2 


/i 2 = nl tan (3' -cos 2/3'., 


(9) 

( 10 ) 


where Mz> = g BB v' 2 (no mixing between U( 1) and £/(1 )b-l is assumed here). From 
-pj), one gets 


sin 2/3' = 


2/r 2 


m 


2 > 
A' 


where m 2 A , = / 1 2 + /r 2 . Note that, again from (9)—(10), we also get 


u' 2 = 


Qj - fA ) - (r 1 + /f|) cos 2/3' 

2g 2 BL cos 2/3' 


( 11 ) 


( 12 ) 


Now we complete our analysis of symmetry breaking. We have 


Vu(v[,v' 2 ) 

= 2/i 2 — 2g 2 BL (v 2 — 3uf), 

(13) 

Vl 2 « v’ 2 ) 

= —2n\ — Ag 2 BL v x v 2 , 

(14) 

V22(v' 1 ,V 2 ) 

= 2gl + 2g 2 BL (3v' 2 -v' 2 ), 

(15) 


where Vq = - g^ X g^ X2 ^ ■ To show that the symmetry will be broken spontaneously, we must 
ensure that the point {v' x ,v'^) = (0,0) is not a local minimum of the potential V. Since 
(V 11 V 22 ~ V 12 ) (0)0) = (2h-i)(2/i 2 ) — (2gl) 2 and Vn(0, 0) = 2g\ > 0 we should impose a 
condition to make (0,0) a saddle point. This condition is 


lA i4 < / 4 - ( 16 ) 

It is worth noting that it is impossible to simultaneously fulhl both the conditions Q and 
(16) for positive values of g\ and ji 2 . However, one should note that the condition (16) 
is valid at the B — L symmetry breaking scale, where the running of the RGEs, from the 
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Figure 1: The evolution of the B — L scalar masses from the GUT to the TeV scale for 
m 0 = Mi/ 2 = A 0 = 1 TeV and Y S3 ~ 0( 1). 


GUT scale down to the B — L breaking scale, may induce negative squared mass of \ 2 . 
At this scale both conditions (J8]) and (16) are satisfied and symmetry is spontaneously 
broken with stable potential. This can be seen as follows. The RGEs for m 2 x 2 are given 
by 


16vr 2 

167t 2 


dm h 

dt 

dm X2 

dt 


12 g 2 BL Ml L - 6 Y/ (ml, + ml. + m| 2j + A( : ,J, 
[‘AQblM'bli 


(17) 

(18) 


where t = In J and Mbl is the gaugino mass of B ', which is given by Mi/ 2 at 
the GUT scale, M X - In order to solve these equations, we should take into account 
all involved RGEs, given in Ref. ra- Rg. 0 reports the result of the running. In 
plotting this figure, we set the following mSugra inspired conditions at the high scale, 
e.g., m 0 = Mi /2 = A 0 = 100 GeV and an order one ratio Y S3 ~ Mr 3 /v'. As can be seen 
from the plot, rn 2 Xi drops rapidly to a negative mass region whereas m 2 2 remains positive. 
Also in Fig. |lj we plot the scale evolution for the scalar mass as well as for and 
m§ 2 ■ The figure illustrates that they remain positive at the TeV scale. Therefore, the 
B — L breaking via a non-vanishing VEV for right-handed sneutrinos does not occur in 
the present framework. 

Before closing this section, let us emphasise that in the BLSSM with Type I Seesaw 
the B — L symmetry is spontaneously broken by the VEV of scalar singlet X 2 only if 
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the Yukawa coupling Y Ur of the term Yv^r^rXi is assumed to be degenerate, i.e. Y Vr = 
Y 0 diag{l, 1,1}. If one assumes a hierarchical texture, where for instance only Y Vr gives 
important contributions, one finds that right-handed sneutrino may acquire a VEV before 
X '2 and breaks both B — L and A-parity. However, here, within the BLSSM-IS, the 
situation is different. One can show that, even with one non-vanishing Yukawa, yy will 
acquire a VEV before the right-handed sneutrino such that B — L is spontaneously broken 
while A-parity remain conserved [15]. 


4 The BLSSM-IS Spectrum 

We have seen the evolution of different parameters from the GUT to the B — L scale. Once 
the B — L symmetry is broken and so is the EW one too (at the Mw scale), different 
particles with different quantum numbers can mix and acquire new mass eigenstates. 
Here, we will focus on the new particles associated with the B — L symmetry, namely, the 
Z' gauge boson, extra Higgs bosons and right-handed sneutrinos. We shall do so in three 
separate subsections. 

Before proceeding to doing so, we should mention that, in all the analyses below, 
we have used the SARAH na and SPheno mm to build the BLSSM. Furthermore, 
the matrix-element calculation and event generation were derived from MadGraph 5 Ed 
and manipulated with MadAnalysis 5 p?]. Finally, notice that all current experimental 

constraints, from both collider (LEP2, Tevatron and LHC) and flavour (BaBar, Belle and 

LHCb) are taken into account in our numerical scans. Also DM and neutrino mass and 
coupling constraints are enforced throughout. 

4.1 The Z' Gauge Boson in the BLSSM-IS 

The U(l)y and U{1)b-l gauge kinetic mixing can be absorbed in the covariant derivative 
redefinition, where the gauge coupling matrix will be transformed as follows: 

G = ( 9yy 9yb ) => G = ( 91 \ , (19) 

\ 9bY 9 BB ) \ 0 9bl ) 

where 


9i = 

9yy9bb 9yb9by 

( 20 ) 

Jg 2 + 9 2 

'BL 

\/g 2 + g 2 , 

( 21 ) 

9 = 

9yb9bB A 9by9yY 

( 22 ) 

\/g 2 + g 2 

Y ^ BB ^ BY 
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In this basis, one finds 


Ml = \(al + alW, M 2 Z , = g 2 BL v 12 + IjjV. 


(23) 


Furthermore, the mixing angle between Z and Z' is given by 

tan 20 ' = 


9i +92 


~9 2 + n v i) 2 9 2 BL-9 2 2 -9l 

4.2 The Higgs Bosons in the BLSSM-IS 


(24) 


m 


The gauge kinetic term induces mixing at tree level between the H® 2 and X 12 states 
the BLSSM scalar potential. Therefore, the minimisation conditions of this potential at 
tree level lead to the following relations 

1 r 


Bn = - 2gg BL v' 2 cos 2/3' + Am 2 H - 4m? H 

O L 

+ (9i + 9 2 + g 2 )v 2 cos 2/3 tan 2/3, 

Bn' = \ [ - 2 g 2 BL v' 2 cos 2/3' + 2 m 2 xi - 2 m\ 2 
+ 99blv 2 cos 2/3 tan 2/3', 


(25) 


(26) 


where tan/3 = ^ and tan/3' = ^f. Note that, with non-vanishing g, the Bji parameter 


depends on v' and the sign of cos 2(5'. We may have both constructive and destructive 


interference between the first term and other terms in Eq. (25). In general, we find that 
the typical value of Bn is of order TeV. 

To obtain the masses of the physical neutral Higgs bosons, one makes the usual redefi¬ 
nition of the Higgs fields, i.e., H ^ 2 = + 2 + 14 * 1 , 2 ) and x \,2 = 75 + + 

The real parts correspond to the CP-even Higgs bosons and the imaginary parts corre¬ 
spond to the CP-odd Higgs bosons. The squared-mass matrix of the BLSSM CP-odd 
neutral Higgs fields at tree level, in the basis (</>i, </> 2 , 4>'ii <^ 2 ), i s given by 


m A,A' ~ 


(27) 


f B/j.ta.n/3 B/i 0 0 \ 

Bn B/j.cot/3 0 0 

0 0 Bfi! tan /3' Bn' 

\ 0 0 Bg,' Bn'cot (3' J 

It is clear that the MSSM-like CP-odd Higgs A is decoupled from the BLSSM-like one A' 
(at tree level). However, due to the dependence of B M on v 1 , one may find m\ = 4^ ~ 


2 _ 2 Bn' 


m A i = 


sin 2/3' 


0(1 TeV). 


r\_/ 













The squared-mass matrix of the BLSSM CP-even neutral Higgs fields at tree level, in 
the basis (ay, cr 2j a '\ ■, & 2 ), is given by 


M 2 


A /, 2 


hH 


V Mf 


T 

hh' 


ML' \ 


K«- 


(28) 


where M 2 H is the usual MSSM neutral CP-even Higgs mass matrix, which leads to the 
SM-likc Higgs boson with mass, at one loop level, of order 125 GeV and a heavy Higgs 
boson with mass ~ m A ~ 0(1 TeV). In this case, the BLSSM matrix M^, H , is given by 


Mh'H 1 


+ 9bl v ? 


V \m 2 A ,s 2 tv ~ 9 2 BL vW 2 


-\™ 2 A'S2P' ~ g 2 BL V 'l V 2 \ 


m 2 A 'S}> + g 2 BLV2 ) 


(29) 


where c x = cos(x) and s x = sin(x). Therefore, the eigenvalues of this mass matrix are 
given by 


m 


h',H' 


1 r 

2 - 


(m 2 A , + Mgi) =F V ( m \> + Ml,) 2 — 4cos 2 2/1' 


(30) 


If cos 2 2/3' < 1, one hnds that the lightest B — L neutral Higgs state is given by 


my ~ 


'm 2 A ,Ml, cos 2 2/3 


1 

f\ 2 


~ 0(100 GeV). 


(31) 


The 


mixing matrix M 2 lh , is proportional 


l 


M hh ' - 


to and can be written as [IS] 
vWi -v\v' 2 \ 

- v 2 v[ v 2 v' 2 ) 


(32) 


For a gauge coupling g BL ~ \g\ ~ 0(0.5), these off-diagonal terms are about one order of 
magnitude smaller than the diagonal ones. However, they are still crucial for generating 
interaction vertices between the genuine BLSSM Higgs bosons and the MSSM-like Higgs 
states. Note that the mixing gauge coupling constant, g , is a free parameter that can be 
positive or negative [T9] . 

In Fig. [2j we show the masses of the four CP-even Higgs bosons in the BLSSM for 
g B L = 0.4 and g = —0.4 . In this plot we fix the lightest MSSM Higgs boson mass to 
be of order 125 GeV. As can be seen from this figure, as intimated, one of the BLSSM 
Higgs bosons, h', can be the second lightest Higgs boson (~ 137 GeV). Both H and II' 
are instead quite heavy (since both m A and m A ' are of order TeV). 
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Figure 2: The BLSSM-IS CP-even Higgs masses versus rriA> for g BL = 0.4 and g = —0.4. 

This sets the stage for the hypothesis made in Ref. [20] (see also [21j), wherein, 
motivated by a ~ 2.9cr excess recorded by the ATLAS and CMS experiments at the 
LHC around a mass of order ~ 137 GeV in ZZ —)■ 4 1 (CMS) and 77 (both) samples 
[22U23U21U25], is was shown that a double Higgs peak structure can be generated in the 
BLSSM, with CP-even Higgs boson masses at ~ 125 and ~ 137 GeV, a possibility instead 
precluded to the MSSM. 

Before proceeding in this respect, though, two remarks are in order: firstly, if g = 0, the 
coupling of the BLSSM lightest Higgs state, h', with the SM particles will be significantly 
suppressed (< 10 ~ 5 relative to the SM strength), so that, in order to account for possible 
h' signals at the LHC, this parameter ought to be sizable; secondly, in both cases of 
vanishing and non-vanishing g , one may fine-tune the parameters and get a light m^, 
which leads to a MSSM-like CP-even Higgs state, H, with ~ 137 GeV. However, it 
is well known that in the MSSM the coupling HZZ is suppressed with respect to the 
corresponding one of the SM-like Higgs particle by one order of magnitude due to the 
smallness of cos(/3 — a), where sin(/3 — a) ~ 1. In addition, the total decay width of 
H is larger than the total decay width of the SM-like Higgs, h, by at least one order of 
magnitude, because it is proportional to (cos a /cos /3) 2 , which is essentially the square of 
the coupling of H to the bottom quark. Therefore, the MSSM-like heavy Higgs signal 
(pp —>• H —>• ZZ —> 41) has a very suppressed cross section and thus cannot be a candidate 
for light Higgs signals at the LHC. 

In the light of this, we will focus in the next section on the lightest BLSSM CP-even 
Higgs, h ', as a possible candidate for the second Higgs peak seen by CMS in both ZZ —>• 4/ 
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Figure 3: The mixing of h ', r 3i , versus the gauge kinetic mixing, g. 


and 77 , see Refs. [Ml [25]. However, before doing so, we ought to setup appropriately the 
BLSSM parameter space, in order to find such a solution. As mentioned above, the recent 
results from CMS indicate a ~ 2.9<r hint of a second Higgs boson at 137 GeV. Herein, for 
definiteness, we consider my = 136.5 GeV as reference BLSSM point. 

As emphasised above, in the BLSSM, it is quite natural to have two light CP-even 
Higgs bosons, h and h', with mass 125 GeV and ~ 137 GeV, respectively. The CP-even 
neutral Higgs mass matrix in Eq. (28) can be diagonalised by a unitary transformation: 

(33) 


r M 2 r* = cliag {ml,m 2 H ,m 2 h ,,m 2 H ,}. 


The mixing elements T 3 2 and P 3 1 are proportional to g and they identically vanish if 
g = 0, as one can see in Fig. [3] Also, in this limit, Tn and approach sin a and cos a, 
respectively, where a is the usual CP-even Higgs mixing angle in the MSSM. 

The lightest eigenstate h is the SM-like Higgs boson, for which we will fix its mass 
to be exactly 125 GeV. As mentioned, numerical scans of the BLSSM parameter space 
confirm that the h' state can then be the second light Higgs boson with mass of 0(137 
GeV). The other two CP-even states, H and H\ are heavy (of 0(1) TeV). The h! can be 
written in terms of gauge eigenstates as 

h' — F 3 i ai + r 3 2 @2 T T 33 a \ + r 34 02- (34) 

Thus, the couplings of the h' with up- and down-quarks are given in terms of the elements 
of the T mixing matrix by 

h'uu: -i h'dd: (35) 

v sin p v cos p 
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Similarly, one can derive the h! couplings with the W + W and ZZ gauge boson pairs: 


ti W + W~ 
h'ZZ 


1 g 2 M w (r 32 sin /3 + r 3i cos (3 ), 

2 4 .9nr sin 2 9' (n^r 32 + 1^31) 

(w 2 r 32 + uir 3 i) ( g z cos 9’ -gsm.9’) 2 


Since sin#' -C 1, the coupling of the h' with ZZ , gh'zz, will be as follows: 


gh'zz — i gz M z (r 32 sin/3 + r 3 icos/3), 

where g z = \Jg\ + g 2 - 


(36) 


4.3 Sneutrino Masses in the BLSSM-IS 


Now we turn to the sneutrino mass matrix. If we write ul,r and S 2 as vl,r = 71 ( 4>l,r + 
i(Jl,r) and S 2 = = (0s + i&s), then we get the following mass matrix for the CP-odd 

sneutrinos: 


m, 


m,-,i = 


<rL<?L 


( 

v ±v 2 v[Re(Y s T Y v *) m 


ml LaR \v 2 v[ Re(YjY:)\ 


m cr L cr R m cr R cr R 


m 


GRCS 


cr<ts 


m 


<rs&S 


(37) 


where m aL<yR) and rn asas are given in Ref. [TO] and are proportional to v 2 , 

vAq, v' 2 , v'fj!, and v 12 , respectively. The mass matrix for the CP-even sneutrino (m S i?) is 
obtained by changing <j LiR) s —>■ 4>l,r,s ■ The sneutrino mass eigenstates can be obtained 
by diagonalising the mass matrices as 


Upiinl 


iUl = m% 


dia 


U^mUUh, = m3 


(38) 


However, the diagonalisation of these matrices is not an easy task and can only be per¬ 
formed numerically. It turns out that the mass of the lightest CP-odd sneutrino, uj, is 
almost equal to the mass of the lightest CP-even sneutrino, v R . Both u IR are generated 
from the mixing between u R and S 2 . The mass of these lightest sneutrinos can be in the 
range 400-1400 GeV, as shown in Fig. [4] 


5 (S)Neutrino Corrections to the Lightest Higgs Boson Mass 

In this section, we calculate the one-loop radiative corrections due to right-handed (s)neu¬ 
trinos to the mass of the lightest Higgs boson when the latter is SM-likc. We show that 
such effects is in the range 400 — 1400 GeV, thereby giving an absolute upper limit on 
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Figure 4: The lightest sneutrino mass as function of mo for = 1.5 TeV and Ao = 2.5 
TeV, so that the SM-like Higgs boson mass is within experimental limits. 


such a mass around 170 GeV ra- The importance of this result from a phenomenological 
point of view resides in the fact that this enhancement greatly reconciles theory and 
experiment, by alleviating the so-called ‘little hierarchy problem’ of the minimal SUSY 
realisation, whereby the currently measured mass of the SM-like Higgs mass is very near 
its absolute upper limit predicted theoretically, of 130 GeV. 

It is important to note that, unlike the squark sector, where only the third generation 
(stops) has a large Yukawa coupling with the Higgs boson, hence giving the relevant 
correction to the Higgs mass, all three generations of the (s)neutrino sector may lead 
to important effects since the neutrino Yukawa couplings are generally not hierarchical. 
Also, due to the large mixing between the right-handed neutrinos IVj and 5V, [25] . all 
the right-handed sneutrinos Dh are coupled to the Higgs boson H 2 , hence they can give 
significant contribution to the Higgs mass correction. In this respect, it is useful to note 
that the stop effect is due to the running of 12 degrees of freedom (3 colors times 2 charges 
times 2 for left and right stops) in the Higgs mass loop corrections, just like in the case 
of right-handed sneutrinos for which there are also 12 degrees of freedom (3 generations 
times 4 eigenvalues). 

To calculate the (s)neutrino corrections to the lightest Higgs mass, we computed, in 
a previous section, the explicit form of the sneutrino masses, while for the neutrino mass 
expressions, which are well known, we refer the reader to nm- Due to one generation 
of neutrinos and sneutrinos, the one-loop radiative correction to the effective potential is 
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given by the relation 
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The hrst sum runs over the sneutrino mass eigenvalues, while the second sum runs over the 
neutrino masses (with vanishing m Vl ). In case of degenerate diagonal Yukawa couplings, 
one finds that the total A V Vj o is given by three times the value of AV^p for one generation. 
This factor then compensates the colour factor of (s)top contributions. 

Therefore, the genuine B — L correction to the CP-even Higgs mass matrix, due to 
the (s)neutrinos, at the scale Q at which d ^ v ”' v ' > = 0, is given by 
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It follows that (see ms for details) 
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Therefore, the complete one-loop squared-mass matrix of CP-even Higgs bosons will be 
given by M t 2 ree + AM 2 , with 

AM 2 = ( ° „ ° „ | , (44) 

VO 8 2 + S 2 )’ y ’ 

where 5 2 refers to the (s)top contribution presented in Eq. (1) of [TOj and 5 2 is the 
(s)neutrino correction given in Eq. (43). In this case, the lightest Higgs bosons mass 
is given by 


m, = 


M\ + Mf + 8 2 + 5l 
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1-4 


M 2 M\ cos 2 2/3 + (<5 f 2 + 51) (M\ sin 2 (3 + M 2 cos 2 (3) 
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For M a M z and cos 2/3 ~ 1, one hnds that 

ml ~ Mf + Si + Si 
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Figure 5: Lightest Higgs boson mass versus the lightest sneutrino mass. 


If m ~ 0(1) TeV, Y v ~ 0(1) and M N ~ 0(500) GeV, one gets that 5 2 u ~ 0(100 GeV) 2 , 
thus the Higgs mass is of order ^/(90) 2 + O(IOO ) 2 + 0(1OO ) 2 GeV ~ 170 GeV. 

In Fig. [5] we present the Higgs mass, m^, as a function of the sneutrino mass, rap, 
for m 0 G [500,1000] GeV and Y Vi couplings varying from 0.1 to 0.4. As can be seen, in 
our model the Higgs mass can be easily within the experimental limits. In particular, 
we have employed the Higgs mass bound as 123 GeV < rrih < 127 GeV [22 EH], where 
we take into account about 2 GeV uncertainty in the Higgs boson mass due to the the¬ 
oretical uncertainties in the calculation of the minimum of the scalar potential, and the 
experimental uncertainties in m t and a s . Finally, notice that all the data points collected 
satisfy the requirement of radiative EWSB. 

6 LHC Signatures of the BLSSM-IS 
6.1 Search for the BLSSM-IS Z' 

Now we study the signatures of the extra neutral gauge boson Z’ in the BLSSM-IS at 
the CERN machine^] The possibility of a Z' decay into a pair of heavy (inert) neutrinos 
would increase the total decay width of the Z'. Therefore, the Branching Ratio (BR) of 
Z l — > l + l~ (l = e,/i), the prime Z' signal at the LHC, is suppressed with respect to the 
prediction of, e.g., the Sequential Standard Model (SSM), which is usually considered as 
benchmark in experimental searches for a Z'. Fig. [6] shows the BRs of all Z' decays. Note 

2 In fact, we assume here that all SUSY particles (including sneutrinos) are large enough so that the 
Z' cannot decay into these, thereby implying that our analysis can be applied also to the standard B — L 
scenario. We will look at some SUSY effects later on in this review. 
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Figure 6: BRs of the Z' decays in the BLSSM-IS as a function of M z > (note that fermion 
species are summed over), for g bl = 0.5 and g = 0.1. 


that we have assumed that the sfermion are quite heavy so that the Z' decay is dominated 
by SM particles and light inert/sterile neutrinos. According to this plot, the BRs of the 
non-SUSY Z' decays are given by |29j 

J2BR(Z' -)■ ll) ~ 16.1%, 

i 

Y BR{Z' -)■ viPi) ~ 7.8%, 

Vl 

Y BR(Z' -)■ qq) ~ 8.92%, 

q 

Y,BR(Z'^u h u h )~ 33.4%, 

Vh 

YBR{Z' —>■ u s P s ) ~ 32.1%, 

Vs 

Y BR{Z' igu h ) ~ 0.6%, (47) 

Vl,V h 

where l, q, and v^, refer to the charged leptons, the six quarks, and the six heavy neutrinos 
respectively. Whereas v s stands for the three inert neutrinos. In this example we have 
assumed = 2.5 TeV, gsL = 0.5, g = 0.1 and heavy neutrino masses are set at 200, 
430 and 600 GeV, respectively. 

It is worth noting that, in our model, the Z' cross sections (c’s) that were used 
to derive the ATLAS and CMS current mass limit could be simply rescaled by a fac¬ 
tor of (gs-L/9z) 2 x (1 — BR(Z' —> new decay channels)). If gs-L = 9z and BR(Z' —> 
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Tabic 2: Representative pp —> Z' —y ee rates (crx BR) for different Z' masses/couplings 


at the LHC (^/s = 8 TeV) in the BLSSM-IS. 


M z > [GeV] 

C^SSM [fb] 

9b-l — 9z 

ct B -l [fb] (with IS) 

= e/ sin 0 W cos 9 W gB-L = 0.5 

9b-l — 0-8 

1000 

170 

6 

41 

105.7 

1500 

21.7 

0.58 

4.5 

13.2 

2000 

3.4 

0.087 

0.72 

2.3 

2500 

0.8 

0.015 

0.15 

0.58 

3000 

0.21 

0.003 

0.04 

0.19 

3500 

0.06 

6 x 10“ 4 

0.009 

0.06 


new decay channels) = 0, this reproduces the SSM cross sections that were used by AT¬ 
LAS and CMS. Considering the scaling of cross sections, the current Z' mass limits will 
be lowered by a factor of a B _ L {Z' —* ll) / (Jssm^' —> U)- This result is consistent with the 
conclusion of Ref. m- 

If M' z = 1000 GeV were considered, BR [Z' —» l + l~) ~ 14% could be achieved (e.g., 
through onsetting Z' decays into inert/sterile neutrinos), in which case a x BR = 16 fb 
when g B -L = 9z = 0.188 and a x BR = 82 fb when Qb-l = 0.5, while in the SSM the 
BR {Z' —y l + l~) ~ 7.6% giving a x BR = 340 fb for both electron and muon channels. 
In this respect, the experimental limit M z > > 2.5 TeV by ATLAS [31j (2.8 TeV |32j by 
CMS) will be lowered, because of a 0.241(0.035] rescaling of the cross section when, e.g., 
g B —L = 0.5(0.188]. This yields a new limit of 1.9(0.7] TeV (2.2(0.81] TeV). For reference, 
Tab. [2] gives o x BR {Z' —> ee) for the SSM and BLSSM-IS at different g B -L values. 



Figure 7: Feynman diagram for qq —>■ Z l —>■ v h v h —> WWll. 

Detecting a Z' signals for such small Z' masses, of order 1 TeV, obtainable for g B -L ~ 
Qz 1 would only be circumstantial evidence for the BLSSM-IS though, as other Z' models 
may well feature a similar mass spectrum. A truly smoking-gun signature of the BLSSM- 
IS would be to produce a Z' and heavy neutrinos simultaneously. Indeed, it turns out 
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that the dominant production mode for heavy neutrinos at the LHC would be through 
the Drell-Yan (DY) mechanism itself, mediated by the Z'. The mixing between light and 
heavy neutrinos generates new couplings between the heavy neutrinos, the weak gauge 
bosons Z, W and the associate leptons. These couplings are crucial for the decay of 
the heavy neutrinos. The main decay channel is through a W gauge boson, which may 
decay leptonically or hadronically. We sketch this production and decay channel via the 
Feynman diagram given in Fig. [TJ Unfortunately, the very fact that gs-L ought to be 
small to comply with current LHC data which rule out Z' detection in the DY induced 
di-lepton channel for 0( 1 TeV) masses in turn means that such Z' —y Vh^h decays also 
remains unaccessible. One needs a stronger gs-L coupling to access the latter, which 
thus requires higher Z' masses. Hence, for the remainder of this studies, we will adopt a 
BLSSM-IS benchmark wherein M Z i = 2.5 TeV. 

Once such a Z' state is produced and decays into zy,Oi —> WWll , one has to further 
sample WW decays. In case of a multi-lepton final state, one ends up with four leptons 
plus missing energy (4 1 + 2 ui), while in case of a multi-hadronic final state states one ends 
up with four jets plus two leptons (Aj + 2l). In addition, it is also possible to have a mixed 
final state (2j+3Z+z/;). (Notice that one or more neutrinos would appear in the detector as 
missing transverse energy/momentum, E™ 1SS -) If two flavours of the heavy neutrinos are 
assumed to be degenerate in mass, one gets the same final states for the produced heavy 
neutrino pair with similar event rates. This will double the number of final state events 
but will make it difficult to distinguish between final state leptons. Therefore, throughout 
the current study, we consider non-degenerate heavy neutrino masses also including the 
interference between every two different ffavours. (See Refs. [33], [33, [35] [36; [3?] 38j for 
alternative phenomenological analyses in the case of the standard B — L model.) 


We thus focus on the possibilities of the LHC in accessing Z' decays into heavy neu¬ 
trinos in the BLSSM-IS. In doing so, we have carried out full Monte Carlo (MC) event 
generation using PYTHIA [39] to simulate the initial and final state radiation, fragmen¬ 
tation and hadronisation effects. For detector effects we have used Delphes [30] . 

We consider the following benchmark: M z > = 2.5 TeV, M Ui = M U5 = 250 GeV, 
M V6 = M Ur = 400 GeV and M Us = M Ug = 630 GeV. Of the three decay signatures of the 
WW pair discussed in Ref. ESI. i.e., 4 j, 2 jlu and 2l2v (l = e), only the latter appeared 
promising, hence we only focus here on this case, by highlighting the main result^] In 
doing so, we produce our results at a/s = 14 TeV assuming a variable luminosity, ranging 

3 A publication with full details and additional results is in progress m- 
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MT 41 [GeV] MET [GeV] 


Figure 8: Number of events versus the transverse mass of the ‘41epton+ Et 1ss ’ system(left) 
and the missing transverse energy (right) The expected SM backgrounds are included . 
The luminosity assumed here is 3000 fb _1 . Note that the bin width is 10 GeV. 

from the standard 300 fb _1 to the tenfold increase forseen at the Snper-LHC J42J. The 
selections assumed in our analysis are as follows [29]: a transverse momentum, pt, cut of 
10 GeV and a pseudo-rapidity, 77 , cut of 2 were set on each electron while the separation 
between two electrons, Ru, was enforced to be 0.2. We have assumed no restrictions on 
E™ 1SS , generated for this signature by two neutrinos escaping detection. 

The key advantage of this channel is that it is almost background free. The main SM 
noise comes from WWZ (three gauge boson) production with a(WWZ) ~ 200 fb at 14 
TeV [121 IHj. In Fig. [8] we show the invariant mass of the ‘4 lepton’ system from the 
Z' signal versus the WWZ background and also the transverse mass of the ‘4 lepton + 
£'™ 1SS ’ system, where such a variable is defined as 

M t = \J(^/M 2 (4l) +p 2 T (4l) + |p?? iss |) 2 - (p T (4/) +p T miss ) 2 . (48) 

These figures indicate that the decay channel ‘4 lepton + _E™ 1SS ’ yields a quite clean 
signature and is rather promising for probing both Z' and after by the end of the 
standard luminosity run of the LHC by using only few simple cuts to extract the Signal 
(S) from the Background (. B ). This is clear from the achievable numbers of events left 
after the set of cuts mentioned above, alongside their statistical significance, S/\/B, as a 
function of the luminosity, see Fig. [9] 


One thing should however be noted at this stage, concerning the statistical analysis 
presented in Fig. [9j Herein, we have at times (i.e., depending on the luminosity) calculated 
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Figure 9: The S and B rates (left) and their significances S/\/B (right), before and 
after the mass cuts, as a function of the luminosity. At the top (bottom) the mass cut 
is M(Al) > 1.5 TeV (M? > 2 TeV). For the significance plots, the colour scheme is as 
follows: the black (red) line is before (after) the corresponding mass cut. 


significances using the expression S/VB even for event numbers of 0(1), which may not 
be entirely appropriate, as such an approach normally requires large event samples and 
B > S. In essence, the S/^/B method is not accurate as it is based on a y 2 distribution 
of test statistics which assumes that noise is negligible. Thus, as an alternative approach, 
we have used the ‘frequentest method’ too, based on Ref. 09- With reference to Sect. 
2 therein, using our MC samples, we have calculated the expected significance at /j = 0 
to reject the null, i.e., background only, hypothesis. For simplicity, we have used only 
the statistical error for the nuisance parameters, which can easily be extracted from our 
histograms. Also, we have computed the expected significance in the signal region before 
and after a somewhat looser final mass cut, in order to allow for a quantitatively sounder 
comparison (e.g., in the case of Mt > 1 TeV only, as the pattern emerging from the M(4Z) 
selection is similar Q From the plots in Fig. 10 it is clear that, at low luminosity, the 

4 Note that the concept of signal region is obvious where the null hypothesis has to be rejected in the 
region where the signal is bumped over the background and be accepted in the control region where the 
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Figure 10: Comparison between the significance defined as S/x/B versus the frequentist 
approach described in the text before (left) and after (right) the cut Mt > 1 TeV. 

difference between the frequentist and S/\/B methods is large while at higher luminosity 
they tend to converge, with the latter approach overestimating the significance over the 
former. Clearly, an appropriate merging of the two approaches is eventually required in 
order to assess the feasibility of our studies at the experimental level, depending on the 
actual size of the data samples. In fact, to do any cut optimisation is out of scope of this 
paper, as in real experimental analyses this is actually done using Multi-Variate Analysis 
(MVA) techniques, which essentially assume a set of chosen cuts for a range of values of 
each of these then calculate the efficiency for each combination. Then various numerical 
algorithms are used to decide on the best choice of combination to finally apply in the 
analysis. What we did here was instead to first characterise the signal and background 
based on MC truth and then decide the relevant cuts after plotting both of these, which 
we have done essentially by chosing visually the region where the former exceeds the 
latter. As we are unable to access such MVA tools, which belong to the repositories of 
the collaborations, we cannot credibly improve our analysis any further at this stage. 


Before closing this section, we should however comment on the influence of the SUSY 
spectrum on the scope of the Z' —>■ v^h signal within the BLSSM. Clearly, once the 
assumption made so far (that all sparticles are heavier than the Z') is dismissed, the Z' 
boson can decay via SUSY objects. This will correspond to an increased value of its total 
width T^, which would then reflect onto the event rates of such a signal. In fact, the 
latter will scale as the inverse of T^U In order to quantify this (reduction) effect induced 
by a low mass spectrum within the BLSSM, we have revisited the benchmarks introduced 

sample is background dominated. 
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Figure 11: The ratio R(Z' ) (as defined in the text) versus the Z' mass for the four heaviest 
benchmarks points in Tab. [2j 


in Tab. [2j excluding only the 1 TeV mass point (now ruled out by data, as discussed) and 
thus including the 2.5 TeV benchmark considered so far in our MC analysis. For each of 
these scenarios we have computed the ratio 


_ r ( z ' all BLSSM decays) 
r(Z' —> SM — like decays) 


(49) 


the inverse of which corresponds to the rescaling factor to be applied to our event rates 
for the Z' —> VhVh signal in presence of a low-lying SUSY spectrum in the BLSSM. We 
can see from Fig. 11 that R(Z' ), obtained from the average total Z' width after scanning 
the entire BLSSM parameter space compatible with current experimental and theoretical 
constraints, is typically smaller than 3, so that the Z' —y rates will generally not 

be smaller than a factor 1/3 with respect to the values considered here. In fact, also 
recall that we have not allowed for decays into muons (l = /x^), which would contribute 
a factor of cs 2 towards the signal rates in the ‘4 lepton + U™ 1SS ’ channel. The SM rates 
presented here would of course be unchanged. One should however include intrinsic SUSY 
backgrounds in the full analysis, as is presently being done in ®t We can anticipate 
that the latter would not spoil the feasibility of such a signal as potentially (i.e., modulus 
experimental verification, both kinematical and statistical) established here over a sizable 
portion of the full BLSSM parameter space. 
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Figure 12: The number of events of the processes pp —> Z —> 2 / 7 * —>• 4/ (blue), pp —$■ h —>• 
ZZ Al (red) and pp —>■ h! —>■ ZZ —>■ 4/ (green) versus the invariant mass of the out 
going particles (4-leptons), 777 . 4 ;. 


6.2 Search for an Extra BLSSM-IS Higgs Boson 


The Higgs decay into ZZ —» 4/ is one of the golden channels, with low background, 
to search for Higgs boson(s). The search is performed by looking for resonant peaks in 
the 7774 ; spectrum, i.e., the invariant mass of the 4 1 system. In CMS m, this decay 
channel shows two significant peaks at 125 GeV and around/above 137 GeV. We define 
by a(pp —> h') the total h! production cross section, dominated by gluon-gluon fusion, 
computed for m h > = 136.5 GeV, for definiteness. (See Tab. I in [20] for the BLSSM 
parameters corresponding to this specific benchmark point, which is well compliant with 
current experimental limits.) From Sect. [4j it is then clear that 


ct(pp -» ti) _ / r 32 \ 2 

a(pp —>• h) SM Ysin/3/ 


(50) 


(wherein the label SM identifies the SM Higgs rates computed for a 125 GeV mass), which, 
for 777/j/ ~ 137 GeV, is of order 0(0.1). Also the ratio between BRs can be estimated as 


BR(h' -> ZZ) 
BR(h ->• ZZ) SM 


( 1 , \ F{Mz/mh') 

V ) F(Mz/m h ) SM 


T 3 i sec/? 


r 32 sin/3 + r 3 i cos/3. 


2 

+ 2 F 


Alw 
m h > 


, (51) 


where 


. 3(1 - 8a: 2 + 20a: 4 ) /3a; 2 -l' 

r (x) = -ttw;—jttw:— arccos 


(4a: 2 - 1 )V2 


2a: 3 


* (2 — 13a: 2 + 47a: 4 ) — |(1 — 6a: 2 + 4a: 4 ) log a: 

Z^L z 


(52) 


First we analyse the kinematic search for the BLSSM-IS Higgs boson, h', in the decay 
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Number of events for 19.7 fb 1 at y/s = 8 TeV 

Higgs mass 

Observed 

(CMS) 

Expected 

(BLSSM) 

Background 

Z -)■ 2 / 7 * 

h -x ZZ 

125 GeV 

25 

18.5 

6.6 

- 

136.5 GeV 

29 

10.2 

9.15 

0.8 


Tabic 3: The observed (by CMS) and expected (from the BLSSM) number of events in 
a mass window around m/, = 125 GeV (121 GeV < m 4 / < 131 GeV) and = 136.5 
GeV (131 GeV < m 4 ; < 152 GeV) in the ZZ —» 4/ channel compared to the expected 
(dominant) pp —> Z —> 2 1'y* —)■ 4/ and pp h ZZ —> 4/ backgrounds. 


channel to ZZ —» 4/. In Fig. [T2j we show the invariant mass of the 4-lepton final state 
from pp — y h' — > ZZ —>■ 41 at = 8 TeV, after applying a Pt[M] cut of 5[2.5] GeV 
on each of the four leptons I2 oQ The SM model backgrounds from the Z and 125 GeV 
Higgs boson decays, pp —Y Z —> 2Zy* —> 41 and pp —>■ h —>■ ZZ —> 41, respectively, are 
taken into account, as demonstrated by the first two peaks in the plot (with the same px 
requirement). It is clear that the third peak at m 4 ; ~ 137 GeV, produced by the decay of 
the BLSSM-IS Higgs boson h' into ZZ —> 41, can reasonably well account for the events 
observed by CMS j24] with the 8 TeV data. This is shown in Tab. [3j where the mass 
interval in m 4 ; that we have investigated to extract the h' signal is wide enough to also 
capture another prominent 145 GeV anomaly seen in the same CMS data set. (The bin 
width used in m is 3 GeV so potentially able to separate the two peaks at 137 and 145 
GeV, yet combining the handful of events under the two peaks makes statistical sense 
given, on the one hand, the small mass difference and, on the other hand, the fact that 
for masses so close to the WW threshold one can find intrinsic BLSSM Higgs widths of 
order GeV.) 

Next we turn to the di-photon channel, which provides the greatest sensitivity for 
Higgs boson discovery in the intermediate mass range (i.e., for Higgs masses below 2Miy)(^J 
Like the SM-like Higgs, the h' decays into two photons through a triangle-loop diagram 
dominated by (primarily) W and (in part) top quark exchanges. As shown in Sect. |4j 
the couplings of the h' with top quarks and W gauge bosons are proportional to some 
combinations of and ^ 2 , which may then lead to some suppression or enhancement 
in the partial width Y(h' —> 77 ). In the SM, BR(/i —> 77 ) ~ 2 x 10~ 3 . Similarly, in 

5 Transverse momentum and pseudo-rapidity thresholds are slightly different for electrons and muons, 
but accounting for this subtlety would not change our conclusions. 

6 The effects of light SUSY particles leading to a possible enhancement of the di-photon signal strength 
of the SM-like Higgs boson were studied in |46j. 
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Figure 13: The number of events of the processes pp h 77 (red), pp —> h! —* 77 
(blue) versus the invariant mass of the outgoing particles (di-photons), m 77 . 


Number of events for 19.7 fb 1 at y/s = 8 TeV 

Higgs mass 

Observed (CMS) 

Expected (BLSSM) 

125 GeV 

610 

666 

136.5 GeV 

170 

177 


Table 4: The observed (by CMS) and expected (from the BLSSM) number of events (after 
subtracting background) in a mass window around m./, = 125 GeV (120 GeV < m 77 < 130 
GeV) and = 136.5 GeV (131 GeV < m 77 < 141 GeV) in the 77 channel. 


the BLSSM, we have found that, for our mh> = 136.5 GeV benchmark, the BR of h' in 
photons amounts to 2.15 x 10 -3 , hence within current experimental constraints. 


The distribution of the di-photon invariant mass is presented in Fig. 13 for a centre- 
of-mass energy yjs = 8 TeV [2D]. Again, here, the observed h —> 77 SM-like signal around 
125 GeV is taken as background while the Z —> 77 background can now be ignored m 
As expected, the sensitivity to the h! Higgs boson is severely reduced with respect to the 
presence of the already observed Higgs boson, yet a peak is clearly seen at 136.5 GeV 
and is very compatible with the excess seen by CMS (23j- This is shown in Tab. [4j In 
fact, ATLAS results point in the same direction as well [22; 23], see |43j and references 
therein. For example, in [22], they have a slightly worse resolution, of 2 GeV, and use a 
mixed sample of 7 TeV (with 4.8 fb _1 ) and 8 TeV (20.7 fb _1 ) data. If we integrate five 
bins around 125(137) GeV, see Fig. 3 of [22], we obtain approximately 700(145) events 
above the background. It is worth mentioning that here we consider both the gluon-gluon 
fusion and vector-boson fusion modes for both h and h' production. 

Before closing this section, we should also mention that the b! —)■ 77 enhancement 
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found in the BLSSM-IS may be mirrored in the 7 Z decay channel EZB for which, at present, 
there exists some constraints, albeit not as severe as in the 77 case. We can anticipate 
(see [21]) that the BLSSM-IS regions of parameter space studied here are consistent with 
all available data. 


7 Conclusions 

In summary, in this mini-review, we have introduced the reader to the minimal SUSY 
version of the well established B — L model with an IS mechanism, that we have termed 
as BLSSM-IS. This scenario nicely combines the theoretically appealing features of SUSY 
with key experimental evidence of Beyond the SM (BSM) physics in the form of neutrino 
masses. 

Initially, we have proceeded with the construction of the BLSSM-IS Lagrangian, fol¬ 
lowed by an illustration of how dynamical EWSB naturally occurs via RGE evolution 
starting from an mSugra inspired model configuration at high scales. Then, we have 
described the emerging particle spectrum, by singling out the dynamics in the three 
specifically BLSSM-IS sectors: i.e., the Z', Higgs and (s)neutrino parts. In three separate 
subsections we have in fact derived the relevant masses and couplings. 

As EWSB and B — L breaking both occur close to the SUSY mass scale of order 1 
TeV, the BLSSM-IS also bears interesting phenomenological manifestations at the LHC 
in the three aforementioned sectors. Therefore, we have studied next the hallmark signals 
of this scenario in turn. Firstly, we described Z' production and decay into a variety of 
leptonic and hadronic signatures preceding via heavy neutrino^J all leading to detectable 
signals at Run 2 of the CERN machine. Secondly, we highlighted the striking feature 
of the BLSSM-IS in the Higgs sector, in the form of a possible additional light Higgs 
resonance yielding sizable 77 and ZZ decays which may even explain some anomalies 
around a mass of 137 to 145 GeV present already in the ATLAS and CMS data of Run 
1 of the LHC and which have not yet been ruled out by current Run 2 samples (neither 
77 [50] nor ZZ [51] ones). 

In short, the BLSSM-IS represents a viable realisation of SUSY, compliant with all 
current data and giving distinctive signatures at the LHC which will enable one to disen¬ 
tangle it from alternative BSM scenarios. These include the simultaneous production of a 
Z' and heavy (s)neutrinos. Further, these can possibly be accompanied by a second light 
CP-even Higgs boson signal (precluded to the MSSM), but not by a 750 GeV one [52, 53] . 
as the BLSSM-IS specific CP-even Higgs state is heavier than its MSSM counterpart (see 

' The case of sneutrino mediated channels has been tackled in Ref. [IS]. 
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[Ml for very recent results in this respect). 
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